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Abstract 

Based on an exact solution of the Bardeen-Cooper-Schrieffer type Hamiltonian on a spherical 
surface, we calculate the specific heat for the electron system with pair correlations on a sphere. 
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We find that it is possible to extract from the specific heat a temperature above which many- 
id ' 

body states with broken Cooper pairs get populated. Therefore, we define this temperature as the 
characteristic temperature signalling the onset of a BCS-type pair-correlated state for electrons on 
a spherical surface. Such spherical electron systems are realized in multielectron bubbles in liquid 
helium, for which the above-mentioned characteristic temperature is found to be of the order of 
10-100 mK. Both the specific heat and the critical temperature show a pronounced (4-6%) odd-even 
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\ parity effect that persists even for numbers of electrons as large as 10 6 . 
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I. INTRODUCTION 



The physics of the two-dimensional electron system remains the subject of intense theo- 
retical and experimental scrutiny Spherical two-dimensional electron systems have not been 
as thoroughly investigated as their flat counterparts, despite the fact that the spherical ge- 
ometry is expected to lead to new physics because of its intrinsic topological difference with 
a flat space p|. Moreover, the presence of curvature also influences the coupling of the system 

n 

to magnetic fields, as has been shown for spherical [2| and toroidal electron systems 
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Spherical two-dimensional electron systems appear in various physical systems: promi- 
nent examples are metallic nanoshells [4] and buckyballs|5]. The most idealized realization 
of a spherical two-dimensional electron system is found in a multielectron bubble (MEB) in 
liquid helium. MEBs are (typically micron-sized) cavities in the liquid, containing anywhere 
from a few up to 10 s electrons. They have first been observed as a result of a surface in- 
stability that occurs when a liquid helium surface is being charged with electrons beyond a 
critical surface charge density, p] 

The bare single electron states on the surface of a spherical bubble are angular momentum 
eigenstates and have discrete energies, characterized by the angular momentum L and with 
degeneracy 2L + 1. At low temperature there will be a well defined Fermi surface located 
at the highest occupied state. Small-amplitude shape oscillations, including surface waves, 
can be quantized as spherical ripplons, described by the Hamiltonian 

#ri P i = ^2 ^Lal^a^m, (1) 

L>0,m 

where a\ m (ol iTO ) is the creation (annihilation) operator for a ripplon with angular momen- 
tum L and its ^-projection m. ujl are the ripplonic frequencies derived in for spherical 
MEBs under external pressure. The electrons can interact with these ripplons. The interac- 
tion between electrons and ripplons on the bubble surface is analogous to that of electrons 
and phonons in solids. In particular, it can lead to the formation of ripplonic polarons. 8j In 
a recent work, the present authors have shown that the electron-ripplon interaction leads to 
the formation of strong pairing correlations below a critical temperature. [9j These pairing 
correlations are similar to the Bardeen-Cooper-Schrieffer (BCS) superconducting pair corre- 
lations in metals. This raises the obvious question of how "superconductivity" will manifest 
itself on a micron-sized spherical surface. More specifically, one can ask how the formation of 
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strong pairing correlations can be detected experimentally, and how the critical temperature 
can be determined. 



In their study of nanoscopic superconducting aluminium particles, Black et al. lOf] used 
electron tunneling spectroscopy to detect the superconducting excitation gap. Also in mul- 
tielectron bubbles, pairing correlations lead to a (pseudo)gap in the density of states. [9| In 
the present paper, we investigate observable manifestations of pairing correlations in MEBs. 
The altered density of states will reflect itself in particular in the specific heat of the system. 
Below, we calculate the specific heat of the spherical two-dimensional electron system and 
argue that the onset of a sharp rise in the specific heat can be used to define a characteristic 
temperature for the transition between a normal state and the pair-correlated state. More- 
over, we derive an analytical expression for this characteristic temperature and show that 
the specific heat reveals a pronounced odd-even parity effect, which persists even for large 
(10 6 ) numbers of electrons. 

II. SPECIFIC HEAT OF COOPER PAIRS ON A SPHERE 

The spherical, two dimensional electron system, is described in second quantization with 
the operators c[ ma and Q >miCT that create resp. destroy an electron with spin a =|, [ in 
the angular momentum eigenstate \£,m). In Ref. {9] we have shown that in the case of 
multielectron bubbles, the electron-ripplon interaction leads to a Cooper-type attractive 
interaction between the electrons, resulting in a BCS-like Hamiltonian: 

00 / 1 

e=o \m=-e a 

-G 2j ct,-m',$,m' ,^t,m,^e,-m,l J ■ ( 2 ) 
m,m'=—£ / 

This Hamiltonian can be mapped onto the so-called "reduced BCS Hamiltonian" studied 



in the context of superconducting nanograms. [U|, [12[ In the present case, represents the 



unperturbed, 2(2£ + l)-degenerate angular-momentum energy level 

where m e is the electron mass and R is the radius of the spherical system. The separation 
between different angular momentum levels i of the unperturbed system is typically much 



larger than the energy of the ripplons, and the interacting electron-ripplon system can 
effectively be treated as a sum over subsystems with different £. The typical energy scale of 
the kinetic energy is 6\ = h 2 /(m e R 2 ). For a fOOOO electron bubble at zero external pressure, 
ei = 0.78 niK. The interaction energy scale G for the same bubble lies around G — 10 mK, 
and the ripplon energy scale for this bubble is 10 /xK. In what follows, we will present results 
as a function of the dimensionless coupling constant G/e\. 
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14) The 



The Hamiltonian (T5]) can be solved analytically, using Richardson's method, 
total energy is a sum of the energies for each £ subsystem, 

oo 

E 3 = J2 E nl 9i ,bV ( 4 ) 

e=o 

and is characterized by the set of quantum numbers j = {rie, ge, ^}^=o,i,2,...- The energy levels 
of the subsystem I are given by 

= (2n, + bt)et - G(n t - g e )(2£ -b e + 2-n e - gi). (5) 

The quantum number ne corresponds to the number of electron pairs in the subsystem, 
whereas bg is the number of unpaired electrons and gi is the number of elementary bosonic 



pair-hole excitations [15|] in the system of rig pairs. In the ground state all electrons are 
paired (pe = 0), except one when an odd number of electrons is present in the system. There 
are no pair-hole excitations = 0) in the ground state. The spherical symmetry, resulting 
in a discrete and degenerate single-particle level structure, allows for an exact solution of 



the Hamiltonian ([2]), describing also the fluctuations 16]. At finite temperatures, two types 
of excitations will be present. The broken pair states (bp ^ 0) correspond on a mean- field 
level to quasiparticle states in the BCS approximation. The presence of collective modes (in 
the context of the Richardson solution sometimes called 'gaudinos') correspond to the case 
when gi ^ |l71 ]. 

At strong coupling (G ^ ei/2) the electron pairs redistribute themselves over approxi- 
mately u = 2G/e\ different I subsystems around the Fermi level £ = Lp to minimize the 
energy. [9(] . The above estimate for the number of partially occupied angular momentum 
levels £ does not hold for G > Nei/2, where N is the total number of electrons in the 

n 

system. As implied from Eq. (25) of Ref. |9|, at G > Ne\/2, when the interaction strength G 
exceeds the Fermi energy for non-interacting electrons, the ground state would correspond 
to an electron configuration with just a single pair per £-level. However, such an extreme 



strong-coupling limit appears irrelevant for the multielectron bubbles under consideration 
and it is not addressed in the following, where the inequality G <C Nei/2 is assumed to be 
satisfied. 

The degeneracy of energy level Ej is a product of the degeneracies of the energy levels of 
the constituent subsystems: 

oo 

h = U J ZeA> ( 6 ) 



£=0 



with 



jit) = 2 b eC 2e+i I 1 9t-0 

where C\ are binomial coefficients. [3] 

The analytic expression for the energy levels and their degeneracies allow to compute the 
specific heat straightforwardly 

where (E) is the statistical average energy of the electrons in the MEB: 

J2j JjEj exp [-Ej/ (k B T)} 
{ } EjJjexvi-Ej/ikBT)} ■ [ } 

The density of states D$(E) allows to express this as an energy integral. The subscript 5 is 
used to indicate that we numerically calculate the density of states by counting the states 
in an energy interval S around the energy E. For the results presented in Fig. [TJ we chose 
energy intervals of 5 = 0.01 e\. Fig. [1] shows the density of states D$(E) as a function of the 
energy above the ground state energy E gs , for different values of the coupling constant Gjt\. 
Each graph shows results for both an even (full circles) and an odd (crosses) number Np of 
electrons on the highest angular momentum level Lp occupied at G = 0. As the coupling 
constant is increased, a clear odd-even effect appears in the density of states. The presence 
of an unpaired electron sensitively raises the degeneracy of the energy levels and thus the 
number of available states. 

The aim of the next section is to examine the effect of pairing interactions, described by 
the Hamiltonian (j2J) and characterized by the coupling strength G, on the behavior of the 
specific heat function of temperature. 
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FIG. 1: The calculated density of states for a spherical 2D electron system with pairing correlations 
is shown as a function of the energy above the ground state energy E gs . Here D$(E) is evaluated 
by counting the number of states in an interval 5 = O.Olei around energy E. Results are shown for 
even (Np = 54, solid circles) and odd (Np = 53, crosses) numbers of electrons on the (unperturbed) 
Fermi angular momentum level Lp = 26. The different panels correspond to different strengths of 
the attractive electron-electron term in the BCS-type Hamiltonian ([2]). 

III. RESULTS AND DISCUSSION 
A. Small coupling, G/e\ < 1/2 

When there are no pairing interactions (G = 0), the energy level spectrum consists of the 
discrete levels q and the specific heat will be similar to that of a collection of fermionic 
quantum rotors. It will be small for temperatures T <C h 2 /(2mR 2 kB) and rise rather 
abruptly to saturate at Nks for T ^> ft 2 /(2mR 2 kB)- The onset of the increase of the 
specific heat at G = is shown in the full curve in Fig. [2j The increase in the specific heat 
starts at a temperature where there appears a non-negligible occupation probability for the 
lowest excited states. These excited states are separated from the ground state by an energy 
~ ei(Lp + 1) and correspond to transitions of electrons between different £ subsystems. 
These transitions will be referred to as interlevel transitions. 
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FIG. 2: The calculated specific heat is plotted as a function of temperature for an even and odd 
number of electrons. The results are shown for several values of the interaction strength G sj ei/2. 
The down and up arrows are positioned at the corresponding characteristic temperatures given by 
(jlip and (|12p for even and odd MEBs, respectively. 

For small, non-zero G <C ei, an additional peak in C(T) develops as shown in Fig. [3 
This additional peak appears at a temperature smaller than the temperature at which the 
G = specific heat starts to increase. This peak corresponds to intralevel excitations of 
the system: pair breaking and pair-hole excitations. As follows from Eq. (J5J), the first pair- 
breaking energy equals the lowest pair-hole excitation energy. For a subsystem £ with an 
even number of electrons, we find Ag = (2£ + 1)G. For an odd number of electrons, this 
becomes = 2£G. 

So, when G <C 6i, the intralevel excitations occur at energies much smaller than the 
interlevel excitations. As G is increased, the specific heat peak that corresponds to these 
intralevel excitations shifts to higher temperatures. In the case of a closed-shell configura- 
tion, no intralevel excitations from the ground state are possible, so that no additional peak 
of C(T) appears at G -C t\. 

B. Large coupling, G/ei > 1/2 

For G > ei/2, the pair-breaking energy becomes larger than the energy spacing between 
the single-electron £ levels. The density of states no longer resembles a set of peaks, but it 
becomes more similar to a step function, where the steps occur at the pair-breaking energy. 
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FIG. 3: The calculated specific heat is plotted as a function of temperature for an even and odd 
number of electrons. The results are shown for several values of the interaction strength G > ex/2. 
The down and up arrows are positioned at the corresponding characteristic temperatures given by 
(jlip and (112p for even and odd MEBs, respectively. 

When the temperature is small enough so as not to populate states with a broken pair, the 
available density of states is small. But when the temperature is increased and broken pair 
states become populated, the relevant density of states jumps to a value that is orders of 
magnitude higher (cf. Figs. QJ). This reflects itself in the behavior of the specific heat: 
as soon as the temperature is large enough to break pairs, the specific heat will be strongly 
enhanced. This abrupt transition can be seen in Fig. [3j 

The values of the temperature that correspond to the fast increase in C(T) clearly cor- 
relate with G. Note that - as distinct from the case of G/t\ < 1/2 - at large coupling 
(G/ex > 1/2) the specific heat is non-zero even at temperatures below the temperature of 
the onset of the rise in C. This is due to the presence of excitations related to interlevel 
transitions of pairs: for large G, these can have energies smaller than the pair-breaking en- 
ergy. With increasing G, the number of these excitations increases (cp. Fig. [TJ: to Fig. [Tfe), 
resulting in the corresponding increase of C with G in the low-temperature region (see 
Fig. 3). 
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C. The critical temperature for pair-breaking 



As implied form the previous discussion and from Figs. [2j [3j both for the case of small 
(G < ex/2) and large (G ^ ei/2) coupling the specific heat starts rising at a temperature 
when states with broken pairs get populated. Below this temperature, BCS-type pairing 
correlations are dominant, but above this temperature, the "superconducting" pairing cor- 
relations are suppressed. So, it is possible to extract from C(T) a typical temperature T* 
separating the two aforementioned regimes. 

Here we make an analytical estimate for T*. The probability for the system to be in the 
ground state is proportional to P gs = J gs exp (— E gs /k B T), where J gs is the degeneracy of 
the ground state. The probability for the system to break up one pair is proportional to 
Pibp = </i& P exp [— (E gs + Al f ) /ksT], since Ap F is the first pair-breaking energy. Here, Jib p 
is the total number of states with one broken pair. An appreciable contribution of states 
with one broken pair to the specific heat appears when P^p becomes comparable to P gs , i.e. 
at a temperature 

kBT ' = ROT (10) 

The ground state for an even bubble is characterized by ge = 0, be = 0. A state with one 
broken pair is characterized by rie —>■ ne — 1 and be — > be + 2. At small coupling, the first 
pair-breaking takes place for £ = Lp. At strong coupling (ei/2 < G <C iVei/2), the electron 
pairs partially occupy about fi = 2G/e\ single-particle i levels around the Fermi level Lp 
and the first pair-breaking can take place on any of these levels (so we still have i « Lp for 
Lp > 1). 

From (JTj) we see that the degeneracy of states that correspond to the first intralevel 
pair-breaking exceeds that of the ground state by a factor J\b p / J gs = 4C2 Lf+1 (~ 8L F for 
Lp ^> 1). The ground state for an odd bubble is characterized by ge — 0, — 1 for the 
lowest partially occupied subsystem with £ m L F , so that the odd bubble has a ground state 
that is a factor 2Cf Lf+1 more degenerate than that of the even bubble. Applying again Eq. 
© we get for odd bubbles Ji bp /J gs = ICf^jCf^ w 8L|/3. 

Given that the pair-breaking energy is equal to the pair-hole excitation energy, why did 
we assume that the increase in the specific heat is due to the breaking of pairs rather than to 
pair-hole excitations ? This is justified because the increase in degeneracy due to the pair- 
breaking is much larger than the increase in degeneracy due to a pair-hole excitation. A 
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pair-hole excitation is characterized by g? — > ge + 1 (whereas bt and ng are unchanged). From 
(J7|) we see that indeed the increase in the degeneracy of states due to the lowest pair-hole 
excitation (for an even bubble) is characterized by the factor (C^ Lf+1 — 1), much smaller 
than Jibpl J gs for Lp 1. 

Thus far, we have only considered intralevel pair breaking. But for G > E\/2, also 
interlevel pair breaking, where an electron pair is broken up and one electron is transferred 
to another i subsystem, become important. 

At strong coupling, the electron pairs partially occupy about \i = 2G/e\ single-particle 
I levels around the Fermi level Lp. All these fi single-particle levels become relevant in 
the pair-breaking process, raising the degeneracy of states due to pair-breaking so that 
Jibp/Jgs ~ ^C^ 2Lp+1 ' ~ 8L 2 F jj 2 for even bubbles and similarly Jibp/J gs ~ 8L 2 f /j, 2 /3 for odd 
bubbles. 



D. Parity effect 

From the preceding discussion, it is clear that there appears a parity effect. The presence 
of an unpaired electron in the odd bubbles alters the degeneracy ratio of the broken-pair 
state over the ground state. The presence of the unpaired electron also changes the pair- 
breaking energy to A Lp = 2L F G as compared with A Lf = {2L F + 1)G for the even bubble. 
The estimate for the temperature T* at which the broken-pair states start being populated 
will therefore also be different for even and odd cases. In particular, we have for the even 
case 



(2Lp + 1) G 
ln(8L|// 

with fi = max[l, 2G/ei], and for the odd case 



y^p -r x;^ 



h T* ~ 2LpG (12) 

kBKdd ln(8LV/3)- {U) 

The location of these temperatures is indicated by arrows in Figs. [21 [3j The estimates 
(Ilip . (ll2j) are in good agreement with the temperatures at which C(T) increases. Note that 
for even bubbles, the rise of the specific heat seems to show an initial "shoulder", and T* ven 
is a bit smaller than T* dd . This shoulder in the specific heat corresponds to the first pair- 
breaking transition. What makes it different from the following pair-breaking transitions? 
As implied from Fig. [1] depicting the density of states, the relative increase in D$ due to 
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FIG. 4: The characteristic temperature above which states with broken pairs become populated is 
shown as a function of the number of electrons on the sphere in panel (a) for an even number of 
electrons. In panel (b) the relative difference between the characteristic temperatures of even and 
odd systems is shown. Here p is the pressure on a multielectron bubble. 



the first pair-breaking transition is larger than that for further pair-breaking transitions. 
Moreover, the energy separation of states with one broken pair from adjacent states is larger 
than that of states with more than one broken pair. 

From Eqs. ffTTI) and (|12p . the relative difference between the characteristic temperature 
for odd and even bubbles is 



V odd even) I e: 



ln(3) 



1 



ln(8L2,/i 2 



2L F + 1 



1 



ln(3) 



ln(8L|// 2 )J 



(13) 



For MEBs with large numbers of electrons (N >> 1), the angular momentum Lp increases 
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with iV approximately as a/ N/2. With increasing N, also the parameter /i = lGjt\ rises 
(it is about 11.6 for N = 10 3 and 170 for N = 10 6 . Therefore, for sufficiently large N, 
inequalities ln(3) <C lia(8L 2 F fi 2 ) -C 2L F + 1 are satisfied so that Eq. f[T3"j) takes the form 

(T; M - T; ven )/T; ven « ln(3)/ ln(8L|/i 2 ). (14) 

The relative difference {T* dd — T* ven )/T* ven , described by Eq. dHj), reflects the fact that - as 
discussed in the previous subsection - for odd bubbles the ratio J\b p / J gs is approximately 3 
times smaller as compared to that in even bubbles. With increasing N, the ratio J\hp/ Jgs 
strongly increases, while the relative difference in Jup/Jgs between odd and even grains 
remains the same at any Lp 3> 1. Since the characteristic temperature T* depends on 
Jibp/ J g s logarithmically [see Eq. fflQ]) ]. the effect of the aforementioned difference in Jib p /J gs 
on (T* dd — T* ven )/T even weakens with increasing N. As a result, the relative difference between 
the characteristic temperature for odd and even bubbles decreases with increasing N, but 
slowly. The characteristic temperature for multielectron bubbles with an even number of 
electrons, is shown in panel (a) of Fig. H] as a function of the number of electrons in the 
bubble and the pressure exerted on the bubble. The pressure on a multielectron bubble 
increases the electron-ripplon coupling, and thus G, and T*. The relative difference between 
the even and odd cases, as given by Eq. f[T3"j) . is shown in panel (b). Even for fairly large 
numbers of electrons, N ~ 10 6 , the T* dd is 4% higher than T* ven . 



IV. CONCLUSIONS 



In this paper, we have focused our discussion on the spherical electron system as it is real- 
ized in multielectron bubbles, even though the results presented here are more general, and 
valid for all the systems well described by the Hamiltonian ([2]). This BCS-like Hamiltonian 
on the sphere can be solved analytically using Richardson's method, and the specific heat 
is calculated from the energy level spectrum and the level degeneracies, or from the density 
of states. The ground state of the system is a state with strong BCS-like pair correlations. 
We show that the specific heat shows a sharp increase as soon as the temperature is large 
enough to populate the many-body states that contain broken pairs. The temperature at 
which the specific heat starts to increase constitutes a characteristic temperature T* that 
separates the BCS-like pair-correlated state below T* from the broken-pair states above T*. 
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This definition relates T* to the critical temperature for the onset of superconductivity that 
can be derived from measuring the specific heat of a superconducting sample. It is impor- 
tant to note that the procedure followed here is not able to distinguish whether the BCS-like 
pair-correlated state appears due to a phase transition or due to a smooth evolution from 
the 'normal' Fermi liquid. The presence of an unpaired electron in spherical systems with an 
odd number of electrons significantly affects both the specific heat and T*: we find that for 
multielectron bubbles with odd numbers of electrons T* is roughly 4-6% larger than in the 
even case. This persistent parity effect is present even for a large total number of electrons. 
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